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AN IMPROVEMENT TO AN ALGORITHM OF 
BELABAS, DIAZ Y DIAZ AND FRIEDMAN 

LOIC GRENIE AND GIUSEPPE MOLTENI 


Abstract. In IBDF081 Belabas, Diaz y Diaz and Friedman show a way to determine, assum¬ 
ing the Generalized Riemann Hypothesis, a set of prime ideals that generate the class group 
of a number field. Their method is efficient because it produces a set of ideals that is smaller 
than earlier proved results. Here we show how to use their main result to algorithmically 
produce a bound that is lower than the one they prove. 


1. Introduction 

We refer the reader to the paper |BDF08| for an outline of Buchmann’s algorithm. 

Let K be a number field of degree tik, with ?t (resp. rf) real (resp. pair of complex) 
embeddings. We denote Ak the absolute value of its discriminant. 

Definition 1. Let W be the set of functions F: [0,+oo) — >• R such that 

• F is continuous; 

• 3e > 0 such that the function F(x)e^ +£ > X is integrable and of bounded variation; 

• .F(O) > 0; 

• (.F(O) — F(x))/x is of bounded variation. 

Let then, for T > 1, W(T) be the subset ofW such that 

• F has support in [0, log T]; 

• the Fourier cosine transform of F is non-negative. 


The main result of IBDFOHj is, up to a minor reformulation: 

Theorem 2 (Belabas, Diaz y Diaz, Friedman). Let K be a number field satisfying the 
Riemann Hypothesis for all L-functions attached to non-trivial characters of its ideal class 
group CL k, and suppose there exists, for some T > 1, an F € W(T) with F( 0) = 1 and such 
that 
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Then the ideal class group of K is generated by the prime ideals of K having norm less than 

T. 


The authors apply the result to the function where L = logT, * is the convolution 

operator and Cl is the characteristic function of (—■§, 4), to get the 
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Corollary 4 (Belabas, Diaz y Diaz, Friedman). Suppose K is a number field satisfying 
the Riemann Hypothesis for all L-functions attached to non-trivial characters of its ideal class 
group CUk , and for some T > 1 we have 
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(Here C = F) k (2k + l) -2 = 0.915965 • • • is Catalan’s constant.) 

Then the ideal class group of K is generated by the prime ideals of K having norm less 
than T. 


Our aim is to find a good T for the number field K as fast as possible exploiting the 
bilinearity of the convolution product. 


2 . Setup 

We use the following definition to simplify a little bit the language. 

Definition 6. A bound for K is an L = log T with T as in Theorem 0 

2.1. Rewriting the theorem. We begin by homogenizing Equation Q and relaxing the 
requirement F( 0) = 1 to F(0) > 0 so that now the condition on the function is 


+oo 


(7) 2 Y log Np Y F<K ^n/ 2 P ^ > F (°) ( lo S Ak - nK -T - n K log(8vr) - 

p m= 1 ” 


+ r 1 


L 


+°° F( 0) - F{x) 
2cosh(x/2) 


dx + hk 


r+oo 

3 2sinh(x/2) 


dx 


Definition 8. Let S be the real vector space of even and compactly supported step functions 
and, for T > 1, let S(T) be the subspace of S of functions supported in 


lo &T log T 
2 ' 2 


Definition 9. For any integer N 1 and positive real 5 we define the subspace S(N,5) of 
S(e 2NS ) made of functions which are constant Mk € N on [kd, (k + 1)<5). 

The elements of S(N,5 ) are thus step functions with fixed step width 5. If N ^ 1, 5 > 0 
and T = e 2NS we have 


(10a) 
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Mk > 1, 
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If, for some T > 1, E S(T ) and F = satisfies (0 then, according to Theorem [2j 

Cl k is generated by prime ideals p such that Np < T. This leads us to define the linear form 
Ik. on S * S by 
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and the quadratic form gx on S by gx(3>) = ^k($ * $)• We can at this point give a weaker 
version of Theorem [2] as 


Corollary 11. Let K be a number field satisfying GRH and T > 1. If the restriction of gx 
to S(T ) has a negative eigenvalue then (Ax is generated by prime ideals p such that Np < T. 

Note that gx is a continuous function as a function from ( S(T ), ||.||i) to R. Therefore if 
logT is a bound for K then there exists an L' < logT such that L' is a bound for K. Note 
also that, in terms of T, only the norms of prime ideals are relevant, which means that we do 
not need the smallest possible T to get the best result. 


Remark. If T > 1 and T E <S(T), then for any e > 0 there exists N 1, 6 > 0 and 
$5 E S(N,5) such that ||<P * $ — ^ e and e 2N5 ^ T. Hence we do not loose 

anything in terms of bounds if we consider only the subspaces of the form S(N, 5). 


2.2. Computing the integrals. Let T > 1 be a real, L = logT and Fj J = Cl* Cl where, as 
above, Cl is the characteristic function of [—4] • We readily see that Fl(x ) = {L—x)C 2 l{x) 
for any x 0. We easily compute 
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where C is Catalan’s constant and dilog(x) is the dilogarithm function normalized to be the 
primitive of suc h that dilog(O) = 0 (this is the normalization of [PARll5| ). 


2.3. A remark on the restriction of quadratic forms. Let q be a quadratic form on 
an ?r-dimensional vector space V of signature ( z,p,m ). We can interpret p (resp. m ) as the 
dimension of a maximal subspace on which q is positive (resp. negative) definite while the 
kernel of q has dimension z = n — p — m. 

Let H be an hyperplane of V and q' the restriction of q to H. A maximal subspace on 
which q' is definite is a subspace on which q is definite, thus the intersection of a maximal 
subspace on which q is definite with H. This means the signature {z 1 ,p' , m') of q' will be such 
that p' ^ p ^ p' + 1 and ml ^ m ^ m' + 1. Cases p = p' + 1, m = m' + 1 and p = p' , rn = rr\! 
are both possible with z = n — p — m = z' — 1 and z = z' + 1 respectively. 
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3. Improving the result 

3.1. Basic bound. We restate |BDF08l Section 3, p. 1191] which determines an optimal 
bound for Corollary 01 Let GRHcheck(K, log T) be the function that returns the right hand 
side of 0 minus its left hand side and BDyDF(K) be the function which computes the optimal 
bound, by dichotomy for instance. The computation of BDyDF(K) is very fast because the 
only arithmetic information we need on K ~ Q [ x \/(P ) is the splitting information for primes 
p < T and is determined easily for nearly all p. Indeed if p does not divide the index of 
Z N/(P) in Ok, then the splitting of p in K is determined by the factorization of P mod p. 
We can also store such splitting information for all p that we consider and do not recompute 
it each time we test whether a given bound logT is sufficient. 

3.2. Improving the bound. We fix a number held K. We denote qn,N,6 the restriction of 
gx to S(N,5). According to Corollary [Til if qK,N,s Las a negative eigenvalue then 2N5 is a 
bound for K. This justifies the following definition. 

Definition 12. The pair ( N,S ) is K-good when gK,iv,<5 has a negative eigenvaliLe. 

We can reinterpret Functions GRHcheck and BDyDF saying that if GRHcheck(K, 25) is neg¬ 
ative then (1,(5) is A'-good and that (l, \ log BDyDF(K)) is A'-good. 

As a first step to improve on Corollary 01 given 6 > 0 we look for the smallest N such that 
(IV, (5) is K- good. Looking for such an N can be done fairly easily with this setup. For any 
i ^ 1, let <f>j be the characteristic function of (—i5,i6). Then (3>i)i^i^iv is a basis of S(N,5). 
We have <Fj * <f>j = F 2 ig = (2 id — |a;|)(|a:|); observe also that the function considered in 
Corollary 0] is j^fTj og 'r. We further observe that 

<f>j * ^j ^\i—j\5 • 

This means that the matrix An of Qk,n,s can be computed by computing only the values of 
Ik (TW) for 1 ^ i ^ 2 N and subtracting those values. 

We then stop when the determinant of An is negative or when 2N5 ^ BDyDF(K). This 
does not guarantee that we stop as soon as there is a negative eigenvalue. Indeed, consider 
the following sequence of signatures: 

(0, p, 0) -A (l,p,0) -A (1 ,p, 1) -A (0,p + 1,2) ->• • • • 

We should have stopped when the signature was (1 ,p, 1) however the determinant was zero 
there. Our algorithm will stop as soon as there is an odd number of negative eigenvalues (and 
no zero) or we go above BDyDF(K). Such unfavorable sequence of signatures is however very 
unlikely and can be ignored in practice. 

The corresponding algorithm is presented in Function NDelta. We have added a limit A r max 
for N which is not needed right now but will be used later. In Function NDelta, we need to 
slightly change GRHcheck to returns the difference of both sides of Equation (0 instead of 0. 
Note that (3>j) is a basis adapted to the inclusion (1 10c: I) so that we only need to compute the 
edges of the matrix An at each step. The test det A < 0 in line Q2] can be implemented using 
Cholesky LDL* decomposition which is incremental. 

One way to use this function is to compute T = BDyDF(K) and for some lV max ^ 2, let 
5 = 2 n*^ and N = NDelta(K, 5, lV max ). Using the inclusion (llOdl) . we see that (N,5) is 
A'-good and that N ^ JV max , so that we have improved the bound. 
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3.3. Adaptive steps. Unfortunately Function NDelta is not very efficient mostly for two 
reasons. To explain them and to improve the function we introduce some extra notations. 
For any 8 > 0, let Ng be the minimal N such that ( N , 8 ) is A-good. Observe that Func¬ 
tion NDelta computes Ng, as long as Ng ^ A max and no zero eigenvalue prevents success. 
Obviously, using (HUcD . we see that for any N ^ Ng, ( N, 8 ) is it-good. We have observed 
numerically that the sequence N 8 ]y is roughly decreasing, i.e. for most values of N we have 
N8 n ^(N + 1)8 n+1 . 

For any A ^ 1, let 8 ]y be the infimum of the S 1 s such that (A, 8 ) is A-good. It is not necessar¬ 
ily true that if 8 ^ Sn then (A, 8 ) is A-good, however we have never found a counterexample. 
The function 8 H > 8 Ng is piecewise linear with discontinuities at points where Ng changes; the 
function is increasing in the linear pieces and decreasing at the discontinuities. This means 
that if we take 0 < 82 < <5i but we have Ng 2 > Ng 1 then we may have Ng 2 82 > Ng 1 8 \ so the 
bound we get for 82 is not necessarily as good as the one for di. 

The resolution of Function NDelta is not very good: going from A — 1 to A the bound for the 
norm of the prime ideals is multiplied by e 2S . This is the first reason reducing the efficiency 
of the function. The second one is that if A max is above 20 or so, the number 8 = log n a^ DF K 
has no specific reason to be near 8 n s ; as discussed above, this means that we can get a better 
bound for K by choosing 8 to be just above either Sn s or 8 i + jy s . Both reasons derive from the 
same facts and give a bound for K that can be overestimated by at most 2<5 for the considered 
A = NDelta(K, 8, A max ). 

To improve the result, we can use once again inclusion fllOdl) and determine a good approx¬ 
imation of 8 n for A = 2 n . We determine first by dichotomy a <5o such that (Ao, do) is A'-good 
for some Ao ^ 1 (we use Nq = 8 in our computation). For any k ^ 0, we take Nk+i = 2A& 
and determine by dichotomy a d/t+i such that (A^+i, 8 k+ 1) is A'-good; we already know that 
is an upper bound for 8^1 and we can either use 0 as a lower bound or try to find a lower 
bound not too far from the upper bound because the upper bound is probably not too bad. 
The algorithm is described in Function Bound. It uses a subfunction OptimalT(K, N,T h T h ) 
which returns the smallest integer T € [Tj, T^} such that NDelta(K, logT/(2A), A) > 0. The 
algorithm does not return a bound below those proved in and . 

3.4. Further refinements. To reduce the time used to compute the determinants, we tried 
to use steps of width 4 8 in [—^ log T, ^ log T] and of width 28 in the rest of [—| log T, | log T], 
to halve the dimension of S(N, 8 ). It worked in the sense that we found substantially the 
same T faster. However we decided that the total time of the algorithm is not high enough 
to justify the increase in code complexity. 

4. Examples 

In this section we will denote T(K) the result of Function BDyDF and Ti(K) the result of 
Function Bound. 

4.1. Various fields. We tested the algorithm on several fields. Let first K = Q [ X \/(P) where 

P = x 3 + 559752270111028720x + 55137512477462689. 

The polynomial P has been chosen so that for all primes 2 ^ p ^ 53 there are two prime 
ideals of norms p and p 2 . This ensures that there are lots of small norms of prime ideals. We 
have T(K) = 19162. There are 2148 non-zero prime ideals with norms up to T(K). We found 
that Ti(K) = 11071 and that there are 1343 non-zero prime ideals of norms up to Tj(K). 
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The time used by Function BDyDF was 58ms on our test computer while the time used by our 
algorithm was an additional 36ms. The test was designed in such a way that our algorithm 
used the decomposition information of Function BDyDF, so it saved a little time. 

We tested also the algorithm on the set of 4686 fields of degree 2 to 27 and small discriminant 
coming from a benchmark of [PARI 15] , The mean value of f° r those fields is lower than 
1 

2 - 

For cyclotomic fields, the new algorithm does not give results significantly better than those 
of Belabas, Diaz y Diaz and Friedman. It might be because the discriminant of a cyclotomic 
field is not large enough with respect to its degree. 


4.2. Pure fields. We computed T(K) and 7\(K) for fields of the form Q[^]/(P) with P = 
x n ±p and p is the first prime after 10 a for a certain family of integers n and a. We computed 
the family of for each fixed degree. The graph shows that it is decreasing with the 

discriminant. The graph of (log log Ak) 2 is much more regular and looks to have a 

non-zero limit, see Figure [Tj below. We computed the mean of -^^-(log log Ak ) 2 for each 
fixed degree. The results are summarized below: 


p 

a ^ 

log A k ^ 

mean 

x z — p 

3999 

9212 

13.19 

x 6 + p 

1199 

13818 

13.38 

21 

x — p 

328 

15169 

13.68 


The small discriminants are (obviously) much less sensitive to the new algorithm. We reduced 
the range for each series to have log Ak ^ 500. The results are as follows: 


P 

a ^ 

mean 

x z — p 

218 

12.35 

x 6 + p 

43 

13.66 

x 21 — p 

10 

17.19 


4.3. Biquadratic fields. We repeated the computations above also for biquadratic fields 
Q[v^T, \/p 2 ] where each pi is the first prime after 10 a * for a certain family of integers a{. We 
found that the mean of -p^y(log log Ak) 2 is 13.63 for the 7119 fields computed and 13.88 if 
we restrict the family to the 1537 ones with log Ak ^ 500. 


Final remarks. In [ BDF08 , Th. 4.3] the authors prove that for a fixed degree T | K) 
(log AKlog log Ak ) 2 and conjecture that T(K) ~ ^P(log AKlog log Ak ) 2 while our computa¬ 
tions suggest that 7\(K) has smaller order. We will prove in a subsequent article (GM15I 
that T(K) x (log AKlog log Ak ) 2 and that Tj(K) <C (log Ak) 2 - 
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Figure 1: -p^j^(loglog Ak) 2 for some pure fields; in abscissa log Ak- 


Input: a number field K 
Input: a positive real <5 
Input: a positive integer JV max 

Output: an N ^ lV max such that (. N,S ) is -ff-good or 0 

1 tab <— (2-/V max + l)-dimensional array; 

2 tab[ 0] <r- 0; 

3 A v- JVmay x N max identity matrix; 

4 N i — 0; 

5 while N < N max do 

6 N 4- N+ 1; 

7 ta6[21V - 1] V- (21V - l)GRHcheck(K, (2IV - l)h); 

8 fafr[21V] A- 21VGRHcheck(K, 2N5); 

9 for i <— 1 to N do 

10 .A [IV, i\ fab[lV + *] — ta&[lV — i]; 

11 A[i,N] 4 - A[N,i]-, 

12 end 

13 if det A < 0 then 

14 | return N; 

15 end 

16 end 

17 return 0; 


Function NDelta(K.<i\lV max ) 
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Input: a number field K 

Output: a bound for the norm of a system of generators of C?k 

1 T 0 <- 4 (log A k + log log A K - (7 + log 27r)n K + 1 + (n K + l) 10 ^ 3 ^ k) ) ; 

2 T 0 f- min (T 0 , 4.01 log 2 A K ); 

3 N <- 8; 5 <r- 0.0625; 

4 while NDelta(K, <5, N) = 0 do 

5 | <5^(5 + 0.0625; 

6 end 

7 T h <- OptimalT(K,iV,e 27V l' 5 - 0 ' 0625 ),e 2Ar5 ); 

8 T <- r h + 1; 

9 while Th < T \ \ T > To do 

10 T <-T h ] N i- 2 N] 

11 T h <- 0ptimalT(K, N, l,T h )\ 

12 end 

13 return T ; 

Function Bound (K) 
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